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Abstract 

We show that pseudo - spin symmetry in nuclei could arise from nucleons 
moving in a relativistic mean field which has an attractive scalar and repulsive 
vector potential nearly equal in magnitude. 
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Almost thirty years ago a quasi-degeneracy was observed in heavy nuclei between single- 
nucleon doublets with quantum numbers (n r , £,}=£ + 1/2) and (n r -l, £ + 2, j= £ + 3/2) 
where n r ,£, and j are the single nucleon radial, orbital, and total angular momentum quantum 
numbers, respectively ■ These authors defined a "pseudo" orbital angular momentum £ 
= £ + 1; for example, (n r si/ 2 , (n r — l)d 3 / 2 ) will have £ = 1 , (n r p 3 / 2 , (n r — l)fs/2) will have 
£ = 2, etc. Then these doublets are almost degenerate with respect to "pseudo" spin, s = 
1/2, since j = £ ± s for the two states in the doublet. This symmetry has been used to 
explain a number of phenomena in nuclear structure || including most recently the identical 
rotational bands observed in nuclei ||. Despite this long history of pseudo-spin symmetry 
|5|||, the origin of this symmetry had eluded explanation. Recenty, it was shown that 
relativistic mean field theories predict the correct spin-orbit splitting ||. In this paper we 
identify a possible reason for this; namely, that the symmetry arises from the near equality 
in magnitude of an attractive scalar, -V s , and repulsive vector, V^, relativistic mean fields, 
V s ~ V„, in which the nucleons move. Such a near equality of mean fields follows from 
relativistic field theories with interacting nucleons and mesons || , with nucleons interacting 
via Skyrme-type interactions [[J, and from QCD sum rules [PH[ . 

A nucleon moving in a spherical field has the total angular momentum j, its projection 
on the z-axis, m, and k = —j3(& ■ L + 1) conserved, where (3 is the Dirac matrix [|ll|]. The 
eigenvalues of k are k = ±(j + 1/2); - for aligned spin (si/2,P3/2, etc.) and + for unaligned 
spin (pi/2)^3/2) etc.). Hence, we use the quantum number k since it is sufficient to label 
the orbitals. The Dirac equation for the single - nucleon radial wavefunction (g K , f K ) in 



dimensionless units is given by [11 



[± + ±±±]g K =[2-E- V(r) ]f K (1) 

[£ + —]/«= [E-A(r)]g K , (2) 
ar r 

where r is the radial coordinate in units of length hc/mc 2 , V(r) = (V (r) + V s (r)) / mc 2 , 
A(r) = (V s (r) - V„(r)) / mc 2 , and E is the binding energy (E > 0) of the nucleon in units 



of the nucleoli mass, mc 2 . First we show that, in the limit of equality of the magnitude of 
the vector and scalar potential, A(r) = 0, pseudo spin is exactly conserved. To do this, we 
solve for g K in (Q) and substitute into ([[]), obtaining the second order differential equation 
for f K , 

+ M+i) + p, (r) _ 2 + B)]/ . = , (3 ) 

ax 1 x ax x 2 

where x = y/E r and 

i— k - 1, k > 0;£ = -k, k < 0, (4) 

which agrees with the original definition of the pseudo-orbital angular momentum []l]|§. For 
example, for (n r s\/2, (n r — l)c/ 3 / 2 ), k = —1 and 2, respectively giving £ = 1 in both cases. 
Furthermore, the physical significance of I is revealed; it is the "orbital angular momentum" 
of the lower component of the Dirac wavefunction. 

Equation (|3[) is a Schrodinger equation with an attractive potential V and binding energy 
2 — E which depends only on the pseudo-orbital angular momentum, £, through the pseudo 
- rotational kinetic energy, ; and not on k. Hence the eigenenergies and eigenfunction 
component f K do not depend on k but only on £. Thus the doublets with the same £ but 
different k (k = £ + 1 and k = —£) will be degenerate, producing pseudo - spin symmetry. 

However, in this limit, there will not be any bound Dirac valence states, only Dirac sea 
states, which contradicts reality. Is it possible that we can have bound valence states and 
quasi- degeneracy for a small A(r)? To answer that question we look at two examples. First 
is the spherical Coulomb potential and the second is the spherical potential well. 

The spherical Coulomb potential for arbitrary scalar and vector fields, V SjV (r) = a s>v /r, 
can be solved analytically JTIJ. The valence eigenenergies are given by, 



4(n + A) 2 [1 - y/(l - (a5)/(n + A) 2 ] - 25(a - 5) 
n ' K ~ (a-5) 2 + 4(n + A) 2 ' [) 

where n is the principal quantum number, n = 1, 2,. . . , a = a s + a v , 5 = a s — a v , and 



A = \k\( Jl + ^2 — 1). The allowed values of k are k = ±1, ±2, . . . , ±(n — 1), — n. The 



dependence on k is in A. If the scalar and vector potential are equal, 5 = 0, then the binding 
energy vanishes, E n K = 0, and hence no bound valence states as stated earlier. For 5 small, 
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E n ,K ~ 7Tt( 1 + ~ FT 21 A 2 7 2\" ' +•••)» ( 6 ) 

2n l \K\n z {4n z + cr) 

Thus the pseudo-spin symmetry is broken in third order in 5. We notice that the breaking 
decreases as n increases and, for a given n, the state with the largest |ac| (which means the 
pseudo - spin partner with k > 0) will have the largest binding energy. Thus pseudo - spin 
quasi- degeneracy coexists with an infinite number of bound valence Dirac states. 

However, the Coulomb potential is not realistic for nuclei and, furthermore, the Coulomb 
potential has higher degeneracies than pseudo-spin since the energies depend only on n, and 
not k in the lowest order. For these reasons we turn now to the spherical potential well: 

V S)V {r) = V S)V >0,r<R; V S)V {r) = 0,r>R; (7) 



The soution of the Dirac equation (0,0) is given in terms of spherical Bessel functions for 

r < R, 

eAk 

9n = A j i+i (z), f K = 2 _ E _y k( z ) > r <R, (8) 



and modified spherical Bessel functions for r > R |]TT 



AK 

9k = A k i+ -(y), f K = -jzTe ' r > R ' (9) 
where z = kr, y = Kr, and the wave numbers are given by 

k 2 = (A - E)(2 -E-V)>0, K 2 = E(2 - E) > 0, (10) 

where e = k/\k\ is the pseudo - helicity, since j = £ + e 1/2, and the eigenvalues of e are ±1. 

The two solutions must match at the boundary, r = R, leading to the two conditions 
which determine the eigenvalues for the same I, but different k: 

Z K j i+1 (Z K ) _ Y K (A - E K ) k /+1 (n) 



ie(z K ) e k k|(y re ) 
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K>0, (11) 



Z^UJJQ _ Y K (A - E K ) k^jYJ 

where Z = kR and Y = KR. Since there are two different equations for the states with 
the same i but different k, the eigenenergies of these two different states will be different in 
general. 

Using the recurrence relations between Bessel functions 

Zj i+1 (Z) = (2£+ 1) jj(Z) - Zy^Z)- Yk i+l {Y) = (2£+ 1) k^Y) + Yk^Y), (13) 

we can eliminate k^_ 1 and rewrite these equations as, 

Z K ] i+l {Z K ) _ r K (A - E K ) k i+1 (Y K ) 



k{Z K ) E K k- e (Y K ) 



k>0, (14) 



Z K j l+1 (Z K ) Y K (A - E K ) k l+1 (Y K ) ~ A 

" jp5 " ^M^) (2£ + 1) ^' k<0 ' (15) 

thereby displaying the fact that the equations become identical for A = producing the 
pseudo-spin degeneracy but, as we shall see, no Dirac valence bound states. 

In Fig.l we plot the left - hand side (LHS) of (|l4] , [l~5D as a function of Z. The LHS 
decreases from a value of zero at Z = to negative infinity at , where Z^~ is the nth 
zero of the spherical Bessel function, jg(Z^~) = 0, with Z = corresponding to n = 0. The 
LHS then becomes discontinuous at this point, and for Z > Z^j it decreases from positive 
infinity to zero at zj°~ +1 , and then negative infinity at zj^J and so on. We call the region 
with Z^°\ < Z < Z , - the nth branch. 

n,l — n+X,K 

On the other hand, the right-hand side (RHS) of both ( |i4|) and ( p~5D increases monoton- 
ically, as illustrated in Fig. 2. The eigenvalues are determined by the points of intersections 
in the nth branch, Z n %, giving the valence eigenergies, 
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n r = n — 1, /c > 0, n r — n, k < 0, (16) 
where the radial quantum number will become clear subsequently. The RHS of ([14]) increases 



from zero at Z = (E = A) to infinity at Z = Z max = y/A(2 -V) R (E = 0) ( see © ). 
However the RHS of (|15D is smaller by an amount (21 + 1)^, and increases from — (2£ + 1) 
at Z = (E = A) to a constant at Z = Z m(l3; (22 = 0). This means that for k < 0, there 
will be a bound state for each branch in Fig. 2 as long as Z n>K < Z max . Furthermore the 
upper component for k < 0, has a zero in each of these branches and thus the radial 
quantum number is then n r = n = 0,1, . . .. However, for k > 0, there will not be a bound 
state for n = since -g = 1 at Z = 0, and thus the two curves intersect only at Z = 0, which 
means = and hence there is no bound state for n = ( see (|8|)). However, there will be 
a bound state for all the other branches in Fig. 2 as long as Z n>K < Z max . Furthermore the 
upper component for k > 0, does not have a zero in the n = ^branch but does have 
a zero in the other branches so the radial quantum number is then n r — n — 1 = 0, 1, . . .. 
This means that the orbit with n r = 0, k < does not have a pseudo - spin partner (in 
fact, this orbital is the "intruder" orbital observed in heavy nuclei), but the orbits with 
n r ,K < 0,n r — 1,k > are in the same branch and are thus pseudo - spin partners which 
agrees with experiment. Also we see from Fig. 2 that the RHS for k > intersects the LHS 
at a smaller Z than k < and thus E„ r _i K> Q > E nr K<0 for the same I in agreement with 
experiment. Furthermore, as the RHS for both k increases, they intersect the LHS at points 
in which the LHS has a larger slope and therefore the points of intersection, Z n)K , are closer. 
Thus these pseudo - spin partners become closer in energy as the radial quantum number 
increases. 

These features can be seen in the limit of a large scalar potential, V s < 1. In this limit, 
the RHS of ( |14|, |l5l ) is large. We use the Bessel function identity fll2| : 



m P ti 1 z-z' o >, z+ z<°> 



which, in the nth branch and for the LHS large and positive, can be appoximated as 
-j i+1 (Z)MZ) » 1/(Z - Z®). If, in addition, is large, 

E nr . l!K>0 « E ( °l , (18) 

where we have denoted E^°\, Y^°~ as the values of E, Y for Z = Z^°], and 



E nr -i,K>o - E nr ,K<o ~ /ta(0) - - "' (0) • (19) 



(2^+1) A 

Thus we see that the energy splitting decreases as the binding energy decreases, which 
is consistent with the fact that pure pseudo - spin symmetry occurs when there are no 
bound Dirac valence states, and that the splitting decreases as the radial quantum number 
increases. Furthermore, for states within the same major shell, the splitting decreases as 
the pseudo - orbital angular momentum decreases. 

Hence, we have shown that pseudo-spin quasi-degeneracy in heavy nuclei can be explained 
by the fact that nucleons in a nucleus move in an attractive scalar, -V s , and repulsive vector, 
V„, relativistic mean fields, which are nearly equal in magnitude, V s ~ V v . The energy 
splitting between states with the same pseudo-orbital angular momentum, £, decreases as 
the binding energy decreases and as I decreases. Although such a near equality of mean 
fields has been derived in specific relativistic field theories f§||, this result probably is a 
general feature of any relativistic model which fits nuclear binding energies, and hence very 



likely a general feature independent of any one model [T(J. In ||, it was shown that V s ~ 
V„ for the isoscalar part of the nuclear mean field (the largest part) but not for the isovector 
part, and the isovector potential has a different shape than the isoscalar potential. This 
implies that pseudo - spin symmetry may be enhanced in heavy proton - rich nuclei, with 
iV ~ Z; these nuclei shall be measured in new radioactive beam facilitiies. 

Pseudo - spin symmetry has been observed also in deformed nuclei [|]]; we are investigat- 
ing the deformed Dirac equation as well. Also, the explanation espoused in this paper implies 
a connection between the wavefunctions of the pseudo - spin doublets. This relationship is 
being worked out. 
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FIGURES 

FIG. 1. The LHS of (0, H) plotted versus Z; Z {0 ) is the nth zero of ji(Z {0) ~) = 0. 



FIG. 2. The LHS of (0, (solid line), the RHS of (g§ (dashed line), and the RHS of 
(short dashed line) plotted versus Z for £ = 1,V = 1.7, A = 0.3, R = 33.5, the radius of 208 Pb in 
dimensionless units. 
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